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The dielectric properties of metallic nanoclusters in the presence of an applied electromagnetic 
field are investigated using non-local linear response theory. In the quantum limit we find a non- 
trivial dependence of the induced field and charge distribution on the spatial separation between 
the clusters and on the frequency of the driving field. Using a genetic algorithm, these quantum 
functionalities are exploited to custom-design sub-wavelength lenses with a frequency controlled 
switching capability. 



Introduction: Recently we developed a theory that de- 
scribes the non-local linear dielectric response of nano- 
metal structures to an externally applied electric field 
[l[. Unlike conventional phenomenological classical the- 
ory [2|, we are able to model the transition from clas- 
sical to quantum response as well as the coexistence of 
classical and quantum response in structures of arbitrary 
geometry. This is of some practical importance because 
metallic nanoclusters can now be made sufficiently small 
such that non-local effects due to finite system size and 
cluster shape dominate the spectral response. In partic- 
ular, when the ratio between the smallest characteristic 
length scale and the Fermi wavelength is comparable to or 
smaller than unity, these systems can fail to fully screen 
external driving fields. Also, in the quantum limit one 
needs to take into account discreteness of the excitation 
spectrum as well as the intrinsically strong damping of 
collective modes. Our model captures these single and 
many particle quantum effects. 

In this paper we report on studies in which we 
explore optimal design of nanoscale metallic structures to 
control electromagnetic field intensity on subwavelength 
scales. We are motivated by the non-intuitive nature of 
quantum response and the potential for applications such 
as surface enhanced Raman scattering [3|, Ij] . 

Model: To capture the single-particle and collective 
aspects of light-matter interaction in inhomogeneous 
nanoscale systems one should consider non-local response 
theory in the quantum regime [1]. The Schrodinger 
equation for noninteracting electrons with mass me and 
charge e moving in a potential V{y)^ is given by 



This equation must be solved simultaneously with the 
Poisson equation to determine the local potential due 
to the spatial distribution of the positive background 
charges. Using the jellium approximation, the resulting 
potential is implicitly given by V^y(r) = 4e7rp(r), where 
the density of the positive background charge p(r) satis- 
fies the condition of neutrality, so that J p{r)dr = Nq\, 
where Nq\ is the number of electrons. The induced po- 
tential 0ind is then determined from the self-consistent 
integral equation 
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Here x(I*^ ^'\^) is the non-local density-density response 
function, and (/)tot(r,cj) = (/)ext(r,cc;) + (/)ind(r,cj) is the 
self-consistent total potential. The induced field is found 
via Eind(r, cj) = — V(/)ind(r,C(;). The external field is as- 
sumed to be harmonic, with frequency Cc;, linearly po- 
larized, and with the wavelength much larger than the 
characteristic system's size. The integral equation is 
discretized on a real-space cubic mesh with lattice con- 
stant L. Its natural energy scale Eq is defined by Eq = 
fi^/(2mei^^), and the resulting system of linear equations 
can be solved numerically. The damping constant which 
determines the level broadening is set to 7 = 2 x 10~^ 

Static Case: To gain better understanding and intu- 
ition into the dielectric response of nanoscale clusters in 
the quantum limit, let us first consider systems consisting 
of two identical nanospheres with a total number of elec- 
trons A^ei and separated by an adjustable distance D (top 
view shown in Fig. 1). The nanospheres are placed in 
a static electric field, and the z-direction of the external 




FIG. 1: The direction of the static external field Eext(ct; = 0) 
points along the z axis. Induced charge density (upper row) and 
corresponding induced electric field logiQ(|Eind(a^, 0, 2;)p) (lower 
row) in systems of two conducting spherical clusters with radii 
R = 7L (shown by black dots), where L is the lattice spacing. 
The plots show top views of the three-dimensional systems. The 
relative distance D between the two closest points of the spheres 
is varied from left to right: (a),(b): D=0; (c),(d): D = 4L, (e),(f): 
D = 13L. A^el = 20 for the system, temperature T=0 K and damp- 
ing constant 7 = 2 x 10~^ Eq. The target volume over which the 
field intensity is to be maximized is indicated as a circle of radius 
L at the system center. 



field is aligned along the line connecting the sphere cen- 
ters. Our goal is to maximize the intensity of the induced 
electric field Wint = Jy |E(r)p(ir in a target volume of 
radius Ry = L, centered between the two clusters, by 
varying the cluster separation D. In the regime of large 
electron densities, classical theory predicts that Wint di- 
verges as the spheres approach each other, i.e. I^int -^ 00 
as D ^ 0, and hence the spherical clusters would need 
to be as close as possible to each other to maximize the 
induced field in the target area. As seen in Fig. 1, this 
is no longer true for small carrier concentrations (here 
A^ei = 20), in which case quantum fluctuations strongly 
influence the electromagnetic response. For sufficiently 
small separations (Fig. 1(a) and (b)) the entire system 
responds as a single dipole (with small corrections at the 
interface between the two clusters). The charge density 
distribution depicted in Fig. 1(a) shows charge polar- 
ization (red: positive, blue: negative) along the applied 
field, whereas the corresponding induced electric field in 
Fig. 1(b) remains relatively homogeneous throughout 
the entire system. Remarkably, as shown in Fig. 1(c) 
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FIG. 2: Landscapes of the induced electric field energy log]^Q(Vl/ind) 
in systems of two conducting nanospheres with radii R = 7L as a 
function of the frequency of the external field lj and the cluster 
separation distance D, measured between the two closest points of 
the nanospheres. The direction of the external field points along 
the line connecting the sphere centers, damping constant 7 = 2 x 
10-^Eo. (a) A^ei = 20, (b) A^ei = 100. 



and (d), there are finite optimum separations between 
the spheres which maximize the induced field at the cen- 
ter between them. As will become apparent, the physical 
reason for this is quantum wave functions constrain acces- 
sibility to geometric features. For the parameters chosen 
in this example, I^opt ^ 4L occurs at the threshold sep- 
aration distance beyond which the two clusters cease to 
act as a single dipole. It is evident from Fig. 1(d) that at 
this resonance condition the overall induced field inten- 
sity is highly inhomogeneous and peaks at two orders of 
magnitude larger compared to off-resonance conditions. 
Moreover, there can be further such resonances, e.g. at 
^opt ^ 7L for the present parameters, which maximize 
the induced field intensity in between the nanospheres. 
As observed in Fig. 1(e) and (f), for larger distances 
D one can ultimately treat the spheres as independent 
dipoles, for which the induced field energy Wint scales as 
D^^ with A = — 8, which we have verified numerically [5|. 
Dynamic Screening: Although the characteristic 
Thomas- Fermi screening length is known to increase with 
decreasing carrier concentration, this quantity only de- 
scribes the screening of slowly varying potentials. The 
lower the carrier concentration, the worse the system 
screens rapidly varying potentials. For relatively small 
electron densities the screening length becomes compara- 
ble with the distance between the spheres D. Thus, sen- 
sitivity of the response of the system to carrier concentra- 
tion suggests strong effects of carrier screening. This will 
be most pronounced in the region between the spheres, 
where the potential undergoes significant changes. To il- 
lustrate how the carrier concentration in the nanospheres 
dramatically changes the dynamic dielectric response of 
the system, we show in Fig. 2 plots of Wint as a function 
of the frequency of the external field and relative distance 
D between the spheres. In Fig. 2(a) we consider the case 



of the same low carrier concentration as in Fig. 1. There 
are A^ei = 20 electrons in the system, with a characteris- 
tic Fermi wavelength Ap ~ 3L, which is the same order 
of magnitude as the radii of the spherical clusters. The 
various observed resonances correspond to excitations of 
different geometric modes available in the discrete spec- 
trum of the system (dipole-dipole, quadruple-quadruple, 
etc). For the parameters chosen in Fig. 2(a), the domi- 
nant geometric resonances occur at frequencies less than 
huj = 0.2^0- Interestingly, there is no zero- frequency 
peak at I) ~ 0, which is in stark contrast to the case of 
denser fillings that correspond to the classical limit, e.g. 
as shown in Fig. 2(b). At low fillings, the delocalized 
charge density response results in less efficient screening 
in the region between the clusters, and hence the system 
of two clusters responds as a whole. This significantly re- 
duces the magnitude of the induced charge densities near 
the closest surfaces of the spheres and limits the maxi- 
mum possible value of VFint- Moreover, quantum dis- 
creteness of the energy levels results in a non-monotonic 
dependence of Wint on D^ which in turn leads to the 
non-zero optimal distance I^opt, discussed above. Note, 
that for the parameter set chosen here at a finite fre- 
quency uj = 0.1 Eq the optimal distance is near D = 0, 
i.e. similar to the static response of the classical system. 
In Fig. El^b) we consider the same system parameters 
but at a higher carrier concentration, i.e. Nq\ = 100 
electrons. The corresponding characteristic Fermi wave- 
length is Ap ^ IL^ which allows the dielectric response 
of the system to be much closer to the classical limit. In 
this case many more geometric resonances are observed 
compared to the low-filling regime (Fig. [SJa)). Also, in 
contrast to the quantum limit these resonances depend 
more strongly on changes in separation D and converge 
into a single dominant peak at hco ^ 0.72Eo for D > 7L. 
Also note the large maximum of Wint dit D ^ and cj ?^ 0, 
as it is expected for the static limit in the classical regime. 
Optimal Design (static): The above example illus- 
trates that there can be significant differences between 
the dielectric response in the classical and quantum 
regimes. Let us now explore how the quantum function- 
ality of such structures can in principle be used for the 
design of nanoscale devices. In the following, we pursue 
an optimal design problem of a prototype system with 
multiple adjustable parameters, using a numerical global 
optimization technique based on the genetic algorithm 
[6]. Specifically, we wish to optimize a system contain- 
ing 5 point-like charges with q = +4e each. In order 
to reduce the complexity of the problem the charges are 
placed on a line along the z axis, and we optimize the 
z coordinates of the placed charges. This reduces the 
optimization problem to 5 parameters. A static exter- 
nal electric field is applied along the z axis. In order 
to discretize the numerical search space, the positions of 
the charges are restricted to be on a lattice with lattice 
constant L. We search for optimal configurations of the 



charges that maximize the induced field intensity in a 
target region Vb, located at the center of the system at 
z = l/2Ltot- Here Ltot is the total length of the opti- 
mization box along the z direction. The total number of 
electrons in the system A^ei = 20 is chosen to insure the 
system's response to be in the quantum regime. It usu- 
ally takes about 1100 seconds to perform optimization on 
a 20 node cluster configured with two IGHz processors 
per node. 
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FIG. 3: Induced field intensity \ogiQ{\'E[n(i(x,0, z)\'^) for different 
configurations of scatterers: five point-like charges with q = +4e 
(marked by white dots), each is placed along a line on the z-axis, 
parallel to the external field. A^el = 20 and 7 = 2 x IO'^Eq. (a) All 
charges are placed at the center, (b) a random configuration, (c) 
optimized configuration with target to yield maximum intensity of 
induced field. 

In Fig. 3(a) the intensity of the induced field is shown 
for the case of all the point charges placed at the center 
of the target area, as it would be suggested by classical 
intuition. While the induced field is indeed largest at 
the system center, the overall intensity is relatively small 
compared to optimized configurations. For purposes of 
comparison, in Fig. 3(b) we also show the intensity of the 
induced field for a random configuration of point charges. 
In contrast to Figs. [3](a,b), Fig. [3](c) displays the induced 
field for a numerically optimized configuration of point 
charges. In agreement with the examples in Figs. 1 and 
2(a), the optimal distance between the placed charges is 
finite. The optimization algorithm finds a compromise 
between the distance to the system center and the inter- 
particle distances of the point charges that maximizes the 
induced charge density. This is achieved via maximiza- 
tion of the induced charge localization in the quantum 
system, leading to the most efficient screening near the 
target volume. Thus, we find that using a genetic algo- 
rithm it is possible to create highly efficient optimized 
structures with broken spatial symmetries which func- 
tion as a sub-wavelength lens for electromagnetic radi- 
ation. Note also that the optimal conffguration in Fig. 
3(c) has an inversion symmetry about the system center 
which arises naturally during the optimization procedure. 
Comparing the optimized result with the classical config- 
uration (Fig. [3Ka)) and the random configuration (Fig. 



[3](b)), it is evident that the optimal configuration leads to 
a field intensity in the focal region that is four orders of 
magnitude larger. We have performed further optimiza- 
tions for the positions of 2-8 point-like charges with the 
same target functionality. Significant even-odd effects are 
observed in the optimal arrangements. For even numbers 
of charges, none of the charges should be placed in the 
system center, whereas for odd numbers of charges the 
optimal configuration consists of two symmetrically ar- 
ranged equal groups of charges, and one charge is placed 
in center of the target area. 
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FIG. 4: 



"Frequency splitter": induced field intensities 



log]^o(l^ind(3^5 05^)| ) for optimized configurations with A^el = 20 
electrons and 10 moving scattering centers (marked by white dots). 
Note that two moving centers are overlapping at z = 5L. We use 
7 = 2 X 10 "^£"0, the external field points along the z direction, (a) 
huJi = 0.09£;o, (b) nu;2 = 0.130£;o. 



Optimal Design (dynamic): Finally, let us consider 
the case of time- varying fields, with the goal to design 
a "frequency splitter" in the sub- wavelength limit. In 
this example we again allow the point-like charges to be 
placed along the z axis, and search for optimal spatial 
configurations of the charges which maximize the induced 
field intensity in a target volume centered at z = 2/3Ltot 
for a field frequency cJi, and in a second target volume 
centered at z = l/3Ltot for a second field frequency 
002' Numerical optimization was performed for 10 mov- 
ing positive background charges with q = +2e each and 
Nei = 20 electrons in the system. In Fig. [H we show 
the induced field intensity for optimized configurations of 
charges with two different frequencies (a) hcoi = 0.09£^o 
and (b) J1LO2 = O.lSO^o- The selectivity of this device can 
be quantified by the induced field energies Wi and W2 
in the target volumes 1 and 2 correspondingly, and their 
ratio at two distinct frequencies uji and UJ2' We find that 
for the optimized configuration the ratio W1/W2 ~ 0.09 
at cji, and W2/W1 ^ 0.07 at uj2' 

Conclusions: 

In summary, the dielectric properties of spherical nan- 



oclusters in their quantum regime offer a richness of func- 
tionalities which is absent in the classical case. These 
include a highly non-trivial screening response and de- 
pendence on the frequency of the driving field. Intu- 
ition based on classical field theory, e.g. divergence of 
induced field in the static limit as the distance between 
the spheres decreases, breaks down, and can thus not 
be relied on for the design of new atomic-scaled de- 
vices. In particular, one cannot expect to induce localized 
charge density distributions with a characteristic length 
scale much smaller than the typical Fermi wavelength 
of the system and collective excitation can be dramati- 
cally modified. Moreover, in the quantum regime the de- 
localized induced charge densities can provide increased 
robustness of the optimized quantity, and thus decrease 
the complexity of optimal design. Quantum mechanical 
effects were also found to set boundaries on the maximum 
values of target quantities, i.e. induced field intensity in 
the system. Using genetic search algorithms, we have 
demonstrated that optimal design can lead to field inten- 
sities orders of magnitude larger than "simple" guesses 
based on intuition derived from classical theory. 

For applications such as surface enhanced Raman scat- 
tering, the electric field enhancement in the nm surround- 
ing the molecule depends critically on the local configu- 
ration of both the metallic particle and the molecule. In 
principle, optimal particle shapes can be found that max- 
imize Raman activity for a given molecule attached to 
the metal surface. The approach discussed in this paper 
opens new possibilities for applications. For example, it 
should be feasible to use our theory in combination with 
search algorithms to create robust and reliable highly Ra- 
man active nano- structured surfaces. 
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